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(S| ■ ABSTRACT 

Ph ■ We present a high angular resolution millimeter- wave dust continuum imaging survey 



of circumstellar material associated with the individual components of 23 multiple star 
systems in the Taurus- Auriga young cluster. Combined with previous measurements 
in the literature, these new data permit a comprehensive look at how the millimeter 
;_( \ luminosity (a rough tracer of disk mass) relates to the separation and mass of a stel- 

t^ ' lar companion. Approximately one third (28-37%) of the individual stars in multiple 

systems have detectable millimeter emission, an incidence rate half that for single stars 
(~62%) which does not depend on the number of companions. There is a strong, posi- 
^ ■ tive correlation between the luminosity and projected separation (op) of a stellar pair. 

iy~. . Wide pairs {ap > 300 AU) have a similar luminosity distribution as single stars, medium 

cn , pairs (op f« 30-300 AU) are a factor of 5 fainter, and close pairs [op < 30 AU) are ~5x 

. I fainter yet (aside from a small, but notable population of bright circumbinary disks). 

^— V I In most cases, the emission is dominated by a disk around the primary (or a wide ter- 

CN ' tiary in hierarchical triples), but there is no clear relationship between luminosity and 

stellar mass ratio. A direct comparison of resolved disk sizes with predictions from tidal 
truncation models yields mixed results; some disks are much larger than expected given 
^ . the projected distances of their companions. We suggest that the presence of a stellar 

C^ , companion impacts disk properties at a level comparable to the internal evolution mech- 

anisms that operate in an isolated system, with both the multiple star formation process 
itself and star-disk tidal interactions likely playing important roles in the evolution of 
circumstellar material. From the perspective of the mass content of the disk reservoir, 
we expect that (giant) planet formation is inhibited around the components of close 
pairs or secondaries, but should be as likely as for single stars around the primaries (or 
wide tertiaries in hierarchical triples) in more widely-separated multiple star systems. 

Subject headings: binaries — protoplanetary disks — stars: formation 



^Hubble Fellow. 



Introduction 



Many, if not most, stars a r e born with close 
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(JAbt fc Levvlll976l : lDuquennov fc Mayor 
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20061 ) ■ Depending on their 



orbits, tidal interactions between individual stellar components in these multiple systems can 
dominate the evolution of their natal circumstell ar material and potentially have drastic con- 
sequen ce s for the planet formation p rocess (e.g., iPapaloizou &: Pringld 119771 : iLin &: Papaloizou 



1979a| ]bl: lArtvmowicz &: Lubowl Il994l ). But even in these hazardo us dynamical environ ments 



many young multiples harbor long-lived circumstellar material (see iDuchene et al.l 120071 ) and a 
growing number of their more mature counterparts are being identified as exoplanet hosts (e.g.. 
Patience et al]|2002l : Elaghavan et aDboodlPesidera k Barbierj|2007l : lDoyle et ahlboill ). Given the 
prevalence of stellar multiplicity, an improved empirical understanding of the dynamical interplay 
between the stars and disks in these systems - including effects like tidal truncation, stripping, 
and the orbital evolution of companions - is fundamental for the development of a comprehensive 
model for the formation of stars and planetary systems. Moreover, constraints on these dynamical 
processes in multiple star systems can be used as high mass-ratio touchstones for theoretical work 
on analogous disk-planet interactions, particularly the creation of tidal gaps and subsequent planet 
migration. 

A wealth of theoretical work suggests that the fate of the circumstellar material in a young 
multiple star system is primarily dependent on the separ ation (o) and mass ratio (q) o f the individual 
components, as well as the orbital eccentricity (e) (e.g.. lArtvmowicz &: Lubowll 19941 ). Systems with 
eccentric orbits have an enhanced likelihood of star-disk tidal interactions. For a given orbit, a 
near equal-mass companion (g ~ 1) should have a more destructive impact on disk material than 
a low-mass companion. But in most cases, the effects of {e, q} on the circumstellar material are 
secondary to the orbital separation. Systems with large separations (a ~ hundreds of AU) should 
impart little or no dynamical effects on their circumstellar material, leaving disks around each 
stellar component that are similar to those around single stars. Conversely, individual disks in a 
small-separation (a ~ a few AU or less) system will likely not survive. Instead, these systems can 
host a circum-multiple disk with a dynamically cleared central cavity out to a radius comparable 
to the stellar separation (~2-3a). However, most multip le systems both in th e field and in young 
clusters have intermediate separations (a ~ tens of AU; iMathieu et al.ll2000l ). The disks in these 
systems may suffer the most dramatic effects of star-disk interactions, resulting in their external 
truncation at a fraction of the component separation (~0.2-0.5a), or their complete dispersal. 

Qualitatively, these theoretical predictions find some observational support. Statistical analy- 
ses of warm gas and dust diagnostics (accretion signatures and/or a near-infrared excess) indicate 
that the presence of a companion wit h separation <40 AU may sign i ficant ly hasten disk dispersal 



near the stars (on ~1-10AU scales; ICieza et al. 



200S 



Kraus et al.l |2012| ) . This diminished fre- 



quency of di s k sign atures for "close" multiples was first noted at (sub) millimeter wavelengths by 
Jensen et al.l (Il994l'l. and later confirmed in surveys o f increasing size and sensitiv ity to dust emis- 
sion (jOsterloh &: BeckwithI Il995l : iJensen et al.l Il996bl : [Andrews &: Williamsl l2005l ) . Since the con- 



tinuum emission at such long wavelengths is primarily optically thin (e.g., iBeckwith et al.lll99d ). 
the systematically lower millimeter-wave luminosities from close multiples compared to systems 
with wider separations or single stars were taken as compelling evidence for decreased disk masses 
due to tidal truncation or disruption. However, that evidence is indirect: those observations relied 
on single-dish photometers that do not resolve the individual stellar components nor their disks. 
A quantitative investigation of the theory of star-disk interactions requires observations that can 
address how disk masses and sizes depend on the properties of the stellar system (particularly a and 
q) . With the right combination of angular resolution and mass sensitivity, interferometric measure- 
ments of the optically-thin millimeter continuum emission from the dusty disks in t hese systems 



are uniquely qualified for that task. Aside from a small collection of systems (e.g., iJensen et al. 



1996a 



Akeson et al. 



199S 



Jensen &: Akeson 



2003 



Patience et al.ll2008l ). such data are rare. 



In this article, we present a Submillimeter Array (SMA) survey of the millimeter-wave con- 
tinuum emission from 23 young multiple star systems in the Taurus- Auriga star formation region. 
These data represent the most comprehensive resolved census of cool dust emission from the disks 
that reside in young multiple systems to date. The motivation for the survey sample is introduced 
in §2, and the observations and data calibration are reviewed in §3. A simple modeling analysis of 
these data is conducted in §4, with a focus on retrieving luminosities and sizes from individual disks 
whenever possible. The results of this imaging survey are synthesized with other information in the 
literature in §4 to extract a statistically representative view of circumstellar material in multiple 
star systems. Based on that analysis, we attempt to reconcile the observations with theoretical 
predictions from tidal interaction models in §5. Finally, our key conclusions are summarized in §6. 



2. The Sample 



Multiplicity searches in the Taurus molecular clouds have a lo ng history of success with a 
variety of techniques, ranging from str a ,ightforward direct irnaging (JReipurth &: Zinnecked 119931 : 



White k Ghez 



ing ( e.g. 



1992, 1995 



2001 



Mathieu et al. 



Correia et al.l 120061 : iKraus &: Hillenbrandll2009bl ) and radial velo city monitor- 



19971) t o more specialized metho ds like lunar occu l tations (ISimon et al 



Richichi et al.l Il999l ) , speckle interfer ometrv (I Ghez et al.l Il993l : iLeinert et al.l Il993l ) 



and most recently aperture-mask interferometry (JKraus et al. 



2011 



There is now a reasonably 
complete census of Taurus multiple systems that have angular separations p ~ 0.03-30", i^-band 
contrast ratios of <6 magnitudes (<4 mags for the systems with the smallest sep arations), an d 



primary spectral types betwe en FO and M4. Assuming a mean distance of 145 pc (ILoinard et al. 



2007: Torres et al. 



20071 . I2OO9I ) and a crude estimation of stellar masses (see lKraus et al.ll201ll ). this 
region of multiplicity parameter-space corresponds to projected separations Up « 5-5000 AU, stellar 
mass ratios q ~ 0.1-1.0 (well into the brown dwarf regime), and primary star masses Mp ~ 0.2- 
2M0. There are currently 71 such multiple "systems" known in Taurus, consisting of 111 "pairs" 
of 179 individual stars. For the sake of clarity, we adopt a simple nomenclature in this article such 
that a "system" refers to any group of associated stars and a "pair" is meant as any subset of the 



system that could potentially interact dynamically. A simple binary (2 stars) counts as 1 system 
and 1 pair. In a higher-order hierarchical multiple like UZ Tau, we consider the 4 stars UZ Tau 
Ea, Eb, Wa, and Wb to comprise 1 system of 3 pairs based on their relative projected separations: 
Ea— Eb, Wa— Wb, and Eab— Wab (i.e., this phenomenological scheme implicitly assumes that pairs 
like Ea— Wa are unlikely to interact directly: all such pairs are listed in §5). 

The selection criteria for our resolved millimeter- wave imaging survey were motivated by prac- 
tical observational limitations and consist of two requirements: (1) a composite system flux density 
of >20mJy at 880 ;um, and (2) an angular separation of >0'.'3 for at least one pair in the system. 
The first criterion is a sensitivity restriction that would ensure that our observations could firmly 
detect (3-5 a) two equivalent disks around individual stellar components with the typical expected 
RMS sensitivity of ~2-3mJy beam~^ (se e ^3). Flux density es t imate s for unresolved systems were 
compiled from the single-dish survey of [Andrews fc Williama (120051). If no suitable 880 //m flux 
densit y was available, the 1.3 mm measurements of iBeckwith et al.l (|l99d ) or lOsterloh &: Beckwith 
(JI995I ) were scaled up by a conservative factor of (1.3/ .88)^ ps 2.2 based on the median emission 



ratio at those wavelengths (lAndrews &: Williams 



2005 



2007al ). Note that for the standard opti- 



cally thin and isothermal assumptions for converting luminosity into mass, this sensitivity threshold 
corresponds to ~5M^ of dust (or a total mass of ~1.5Mjup for a 100:1 gas-to-dust mass ratio). 
The second criterion is a resolution restriction set by the longest available baselines of the SMA 
(~0.5 km) that ensures we would be able to resolve the individual stellar components of a given pair. 
Only systems where the angular separations of all its constituent pairs are <0'.'3 were excluded. 
The resulting sample includes systems where some pair has a projected separation ap > 40 AU. 



Of the original 71 systems, 29 exceed the 880 /im luminosity selection threshold. Of those 29 
systems, only 2 fail to meet the resolution criterion (IS Tau and DQ Tau; the latter is a spectroscopic 
binary). The resulting 27 systems are comprised of 52 pairs and a total of 77 individual stars. We 
elected not to observe 4 of those systems with the SMA (FS Tau/Haro 6-5B, XZ/HL Tau, FZ/FY 
Tau, and V807/GII Tau) because their wide-separation pairs were already resolved with single-dish 
telescopes and their remaining pairs were too faint or too close to meet our selection criteria. The 
systems in our 880 fim flux- and resolution- limited sample are listed in Table [H 



Given the practical restrictions that were imposed in its construction, it is important to in- 
vestigate the resulting sample for unintended biases. Since our selection criteria do not specifically 
address the stellar masses in these systems, the key potential biases could be related to the mass 
ratios of pairs (g), the "primary" masses {Mp; meaning the mass of the brighter component of 
a pair), or the total pair masses (the sum of the primary and secondary components of a pair. 



M, 



pair 



Mp + Ms). Figure [U directly compares the cumulative distributions of q, Mp, and Mpair for 



the pairs included in our sample (black) against those that were excluded due to their low millimeter 
luminosi ties and /or small an gular separ ations (gray) . Not e that masses were typically determined 
from the lBaraffe et al.l (|l998l ) models by lKraus et al.l (|201ll ): complete references are provided in §5. 
A two-sided Kolmogorov-Smirnov (K-S) test confirms that our sample is not biased with respect to 
q or Mp, with probabilities that the sample multiples and other multiples are drawn from the same 
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Fig. 1. — A comparison of the cumulative distributions of the stellar mass ratios (q), primary masses 
{Mp), and composite pair masses (Mpair = Mp + Ms) between the systems that were selected {black) 
and excluded (gray) from our sample. While our selection criteria do not produce a significant bias 
in q or Mp, they tend to include pairs with marginally higher total stellar masses (Mpair)- 

parent distributions of mass ratios or primary masses being ~73% and 20%, respectively. However, 
the distribution of Mpair for the pairs in our sample are found to be drawn from a marginally 
different ~ and systematically higher - parent distribution than their counterparts that were not 
selected, with a K-S probability of only 2% (although it is worthwhile to keep in mind that these 
individual stellar mass estimates are crude). 



3. Observations and Data Reduction 

The 23 multiple systems listed in Table [1] (not in cluding the add itional 4 systems in italics; 
see §2) were observed with the SMA interferometer (JHo et al.l |2004| ) in a variety of observing 
configurations and receiver settings over the past 6 years, with most of the data obtained in the 
past 20 months. An observing log is provided in Table EJ All systems were observed in the compact 
(C) array configuration, with baseline lengths of 8-70 m. Additional measurements were made with 
the extended (E: 28-226 m baselines) or very extended (V: 68-509 m baselines) configurations for 
the systems that contain pairs with smaller aiigular separatioiis. Som e of these data were presented 
in previous work ( Andrews fc Williamsll2007bl : [Andrews et alj|201ll ). All but 4 of these 23 systems 
were observed with the 345 GHz receivers: the pairs in the HP Tau, GI/GK Tau, MHO 1/2, and GG 
Tau systems that we aimed to probe have wide enough separations that they were instead observed 
with slightly lower resolution at 230 GHz. We made some additional observations of systems that 
were not in our sample, since they had not yet been observed at millimeter wavelengths (see §5). 

In most cases, the SMA dual-sideband receivers were tuned to a local oscillator (LO) frequency 
of 340.755 GHz (880 ^m) or 225.497 GHz (1.3mm). Some tracks used shifted LO settings to acco- 
modate other projects that shared one night of observing. The data obtained in 2010-2011 employs 
two IF bands (per sideband) spanning ±4-6 GHz and and ±6-8 GHz from the LO frequency (only 
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Fig. 2. — Mosaic of SMA continuum images for wide-separation multiple systems. Contours are 
drawn at 3 o" intervals (see Table ED; hatched ellipses mark the synthesized beam dimensions. Stellar 
positions are marked with blue (primary), red (secondary), or purple (intermediate mass in higher- 
order systems) crosses. A blue dashed curve marks the SMA primary beam in the FV Tau map. 
The FV Tau inset shows the data synthesized at higher resolution, with a beam size of Gf!7 x C/.'S. 
The panels are drawn to scale; a 500 AU scale bar is marked in the HP Tau system panel {top left). 



the lower IF band was available for the 2 observations in 2005 and 2009). Each IF band contains 
24 partially overlapping 108 MHz- wide spectral chunks (per sideband). Aside from one chunk re- 
served for the local CO transition, each of these was coarsely divided into 32 channels to observe 
the continuum. A finer sampling of 256 channels per chunk was used to probe the CO emission, 
corresponding to a velocity resolution of 0.40 and 0.55 km s~^ near the J=3— 2 and J=2— 1 transi- 
tions, respectively. The observations cycled between various target systems and the nearby quasars 
3C 111 and J0510-I-180 on timescales of ~20 minutes for the compact array and 10 minutes for 
the longer baseline configurations. Bright quasars (3C 279 or 3C 454.3), Uranus, and satellites 
(Titan, Callisto) were observed as bandpass and absolute amplitude calibrators when the targets 



were at low elevations. Observing conditions were often excellent, with precipitable water vapor 
levels ranging from 1.0-1.6 mm and ~2-3mm for the 880 //m and 1.3 mm observations, respectively. 

The data were reduced with the IDL-based MIR software package. The spectral response was 
calibrated using observations of bright quasars as references, and the central 82 MHz from the in- 
dividual spectral chunks in each sideband and IF band were averaged into an effective continuum 
channel (excluding the chunk containing a CO transition). The antenna-based complex gain re- 
sponse of the system was determined using the phase calibrator nearest to the target. The absolute 
amplitude scale was set based on observations of Uranus or planetary satellites, and is expected to 
be accurate at the level of ~10%. For each target, the continuum channels for both IF bands and 
sidebands from each set of observations were combined into a composite set of calibrated visibilities. 
The MIRIAD software package was used to Fourier invert those visibilities, perform a deconvolution 
using the CLEAN algorithm, and restore the CLEANed maps with a synthesized beam. The synthe- 
sized beam dimensions and RMS noise levels for the naturally weighted datasets are provided in 
Table El The SMA continuum maps are shown in Figures [2]|11 In most cases, the observations of a 
given multiple system were not sufficient to clearly detect CO emission from any circumstellar gas: 
the few exceptions will be discussed elsewhere. 



4. Disk Properties from Simple Emission Models 

The SMA survey observations described above comprise the largest resolved millimeter-wave 
census of circumstellar material in young multiple star systems to date. In this section, we aim 
to measure two fundamental properties from these data - luminosities (which are related to dust 
masses) and sizes - for the disks around the individual stellar components in each multiple system. 
These basic disk parameters are estimated by fitting a simple model of the continuum emission 
morphology directly to the observed visibilities. The 27 multiple systems in our sample contain 
77 individual stars. The available data have sufficient angular resolution to associate any dust 
emission with 50 of those stars. The individual components in the close pairs MHO 2 AB, T Tau 
Sab, FS Tau AB, DI Tau AB, UX Tau Bab, XZ Tau AB, GG Tau Aab, UZ Tau Eab, V807 Tau 
AB, GH Tau AB, HP Tau/G3 AB, HV Tau AB, and Haro 6-37 Aab are not resolved. We treat the 
millimeter signal from each of these 13 pairs as if it arises from a "composite" disk (see §5.1). 

Including two of those composites, there are 14 individual disks in this sample that are suffi- 
ciently well-resolved to provide robust estimates of their basic parameters. In these cases, we define 
a simple, azimuthally symmetric and geometrically flat emission model with a power-law radial 
surface brightness distribution, I^ oc R~^, that extends to an outer edge, R^- The emission profile 
is normalized such that the total flux density F^i = J Ii^dQ. This parametric emission morphology is 
designed to mimic what would be expected from a disk structure model with power-law surface den- 
sity {Tid ex R~^) and temperature (T^ ex R^'^) profiles. Pressing that resemblance, the radial index 
X is analogous to the sum p+q and the normalization F^ is a rough proxy for the product Kd{Td)M(i, 
where Hd is the dust opacity, (Td) is a characteristic temperature, and Md is the dust mass, mod- 
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Fig. 3. — Same as Fig. [21 for medium separations. 
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Fig. 4. — Same as Figs. [2] and [3l for small separations. 
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ulo small correction factor s for any high optical depths in the disk center (JBeckwith et al.lll99Cl : 
Andrews fc Williamsll2005l ). Our data do not generally have enough sensitivity on long baselines to 
provide useful quantitative constraints on both the emission g r adient and size, which are eff ectively 
degenerate at this modest resolution (see iMundv et al.l 119961 : [Andrews &: Williams! l2007bl ) . Since 
the key parameters of interest from the perspective of tidal interaction models are {F^, Rd}, we 
elect to fix the gradient to a fiducial value, x = 1.5, motivated by the standard assumptions for 
irradiated accretion disks {p = 1, q = 0.5; see lHartmann et al.lll998l ). For reference, adjustments to 
the radial index of ±30% (Ax ~ ±0.5) induce systematic changes in the size estimates of ~20-40%: 
steeper (shallower) gradients produce larger (smaller) sizes. 

In addition to the two free parameters in the surface brightness model, {F^, Rd}, there are 
formally five other parameters related to the projection of the model into the sky plane: the disk 
center relative to the observed phase center {Aa, A6}, the disk viewing geometry described by its 
apparent inclination and orientation {id, PAj^}, and the distance to the observer {d| . The latter is 
fixed to d = 145 pc, with a sy stematic uncertainty estimated to be roughly ±10% (JLoinard et al 



2007 



Torres et al.l 120071 . l2009l ). We fix the centroid positions before estimating other parameters, 



typically based on an elliptical Gaussian fit to the visibilities for individual components that exhibit 
continuum emission. In general, that technique recovers the expected stellar positions well within 
the position accuracy of the SMA data (~0'.'l in an absolute sense, and considerably better in a 
relative sense for the few cases with multiple disk detections). For stellar components that do not 
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Fig. 5. — Model fit demonstration for the HK Tau binary. From left to right, the panels show the 
SMA 880 ^m image, the best-fit model image, and the imaged residuals. Contours are drawn at 
3mJy beam~^ (3o") intervals. The synthesized beam is shown in the lower left of the data panel. 
The HK Tau A and B stellar positions are marked with blue and red crosses, respectively. 



exhibit any millimeter emission or may be partially blended with other components, we rely on 
the positions (or projected angular separations and orientations) provided from optical/infrared 
measurements in the literature to assign their {Aa, A6} values. In practice, we estimate the 
best-fit values of 4 free parameters {F^, Rd, id, PA^} and their uncertainties for each resolved 
disk by comparing model p redictions directly with the SMA visibilities using the non-linear x^ 
minimization routine MPFIT (JMarkwardtll2009l ). In each case, several randomized initial parameter 
sets were employed to avoid trapping in local minima. The results are compiled in Table [H Figure 
[5] shows an example model fit for the well-separated and resolved disks of the HK Tau binary. 

For 4 of the 14 disks with modeling results in Tabled the emission is not sufficiently resolved to 
place meaningful constraints on the disk viewing geometry. In those cases (FX Tau A, UZ Tau Wb, 
HN Tau A, and Haro 6-37 B), we assumed a fiducial {id = 45°, PA^ = 90°} in order to calculate 
a reasonable estimate of {Fd, Rd} (alternative viewing geometry selections produce the same flux 
densities and sizes within the quoted uncertainties). The viewing geometries for the HK Tau B, HV 
Tau C, and GG Tau Aab disks were fixed based on observations of their scattered ligh t morphologies 
or molecular kinematics (jPuchenelboioljMcCabe et al.|[201ll : lGuilloteau et al.|[l999l ). The UX Tau 
A and GG Tau Aab disks were modeled as rings, with empty central r egions out to radii of 25 and 



185 A U based on the more sophisticated analyses of [Andrews et al.l (J201ll ) and iGuilloteau et al. 



(|l999l ). respectively. The UZ Tau Eab circumbinary disk was also modeled as a ring, with no 
emission inside a radius of 15 AU: a detailed analysis of this disk will be provided elsewhere (Harris 
et al., in preparation). Of the remaining 49 individual stars or close-pair composites, 27 have 
firmly detected - but unresolved - millimeter emission and 22 others do not. For simplicity, point 
source models were used to measure Fd for the population of detected, but unresolved, disks. After 



11 



subtracting models of any emission from nearby stars, upper limits (3 a) on F^ were estimated for 
the undetected components by computing the RMS noise level in a 4" x 4" box centered on the 
stellar position in a synthesized residual map. The point source flux densities and upper limits for 
these 49 individual stellar components (or close-pair composites) are compiled in Table [5l 

The emission from the individual disks in the RW Aur and DD Tau binaries is blended. In 
these cases, we adopted an iterative modeling strategy. First, a single disk model was fitted to the 
dominant emission component (the primaries) and subtracted from the data. An initial estimate of 
the disk model for the blended component (the secondaries) was made from the residuals. Based on 
those results, a composite model with both disk components was then used to fit the data and derive 
proper parameter estimates and uncertainties. This method naturally accounts for the blending 
with inflated formal parameter uncertainties for each model component. None of the disks in these 
pairs is spatially resolved, so the underlying model for each component disk is a point source. 



5. Results 

We aim to take advantage of this resolved millimeter-wave census to address some key aspects 
of disk evolution in the presence of stellar companions. To do that, we link our SMA survey 
results into a comprehensive compilation of stellar separations (p), component masses {Mp and 
Ms) and mass ratios (q), and millimeter luminosities (defined as the summed emission in a pair: 
-^pair = Fd,p + Fd,s) for all of the known potentially interacting stellar pairs in Taurus (with spectral 
types FO to M4). Those data are compiled in Table [6l A complementary list of single stars in 
Taurus with available millimeter-wave observations in the literature is provided in Table [71 We 
should note that the re are 36 other single stars in this spectral type range in the compilation of 



Luhman et al.l (120101 ) that, to our knowledge, have not yet been observed at millimeter wavelengths. 



5.1. Millimeter Detection Statistics 

If interactions with companions efficiently remove material from circumstellar disks in young 
multiple systems, there should be a clear signature in the relative detection fractions of millimeter- 
wave emission (a rough proxy for dust mass) between isolated (single) stars and the individual stellar 
components of multiple systems: the fraction of stars that exhibit detectable millimeter emission, 
/mm = A^det/A'tot) should be substantially h igher for singles compared to multiples. This feature 



has been noted anecdotallv in the past (e.g.. ljensen et al.lll994l . Il996bl : lOsterloh fc Beckwith 



1995 



Andrews fc Williamsll2005l ) . but the inability to assign millimeter emission to individual components 



in multiple systems has limited any firm quantitative assessment of the detection statistics. 

There are millimeter-wave continuum measurements available for 52 single stars in Taurus, 
48 binaries (96 stars), 13 triples (39 stars), 7 quadruples (28 stars), 2 quintuples (10 stars), and 1 
sextuple (6 stars; the LkHa 332 system). Using the component-resolved F^ measurements in Tables 
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m and [5] along with the additional literature photometry compiled in Tables [6] and O we computed the 
/mm values for single stars and multiples listed in Table [51 The ranges of A'^^et and /mm for multiples 
correspond to the potential distribution of the observed millimeter emission among any unresolved 
components. We find that 32 of 52 single stars show millimeter emission, /mm = 62 zt 11%, while 
only 50-67 of 179 individual stars in multiple systems have millimeter-wave detections, /mm = 
28-37±5%. A two-tailed Fisher Exact test confirms that these are indeed statistically different 
detection fractions, with a p-value < 0.002. Remarkably, the millimeter detection fraction for 
individual stars in multiple systems does not depend on the number of companions: binaries, 
triples, and higher-order groups all have /mm ~ 1/3, roughly half the detection rate for singles. 
That uniformity is a testament to the hierarchical nature of multiples, where higher-order (N^, > 3) 
systems tend to be be constructed of sets of binary pairs. 

Given that the detection fraction for stars with companions is roughly half that for stars with- 
out them, it may seem natural to assume that only one stellar component in a multiple retains disk 
material. While not uncommon, this is not necessarily the typical scenario. Of the 48 binaries in 
Taurus, 20 exhibit millimeter emission. Of those 20 pairs, our millimeter observations have resolved 
the individual components of 12 (from 24 stars). In 6 of those pairs, the emission is concentrated 
solely around the primary. In the other 6, both components show some emission - however, the 
primary is always brighter. So, for the sub-population of component-resolved binaries with mil- 
limeter emission, the detection fraction for individual stars is actually fairly large (18/24). In the 
higher order multiples, the situation is slightly more complicated by their hierarchical structure. 
In some cases, we find dust emission coincident with all components, albeit usually with some pair 
presumably surrounded by a circum6mar?/ disk (e.g., UZ Tau, MHO 1/2). In others, the dust 
emission is only present around a more isolated, distant companion (e.g., HV Tau, Haro 6-37). 

As might be expected, the likelihood that any individual component of a multiple system 
harbors a circumstellar disk that is massive enough to generate detectable millimeter emission 
depends critically on the individual details of the system. The following sections explore the 
potentially observable signatures expected from star-disk interactions, with a more explicit focus 
on some key connections between the stellar properties tied to orbital dynamics and the basic disk 
characteristics that can be inferred from the millimeter data. 



5.2. Pair Demographics and Disk — Star Connections 

Theoretical models of star-disk interactions in binary pairs suggest that the separation be- 
tween the stellar co mponents is the key property that controls the tid al truncation of individual 



circumstellar disks (jArtvmowicz &: Lubowl Il994l : iPichardo et al.l |2005| ) . These interactions effec 



tively remove mass from the outer regions of these disks (either by accretion or ejection into the 
local interstellar medium), such that stellar pairs with smaller separations should harbor smaller - 
and therefore less massive - disks. Since the cool dust in the outer disk that is being stripped by 
this process emits continuum radiation at millimeter wavelengths, these tidal interactions should 
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Fig. 6. — (a) A comparison of millimeter flux densities from potentially interacting pairs as a 
function of the projected pair separation. Single stars are shown to the right of the plot as 
black points for reference. The pair population can be distinguished into 4 clear sub-categories: 
wide (cp > 300 AU), medium (cp = 30-300 AU), and small (cp < 30 AU) pairs, and circumbi- 
nary disks (purple), (b) The cumulative distributions of millimeter flux densities for each of these 
sub-categories, constructed with the Kaplan-Meier product-limit estimator to include the available 
upper limits. Millimeter luminosity is strongly dependent on projected separation. 



naturally produce an observable trend where the millimeter luminosity is positively co rrelated with 



the separatio n of the stellar pair. Indeed, the pioneering work on this subject by iJensen et al. 



(jl994l . Il996bl ) clearly identified that the millimeter luminosities for pairs with projected separa- 
tions Up < 50-100 AU were statistically lower than their more widely-separated counterparts or 
single st ars. These different popula t ions were confirmed with l arger and deeper millimeter-wave 
surveys (lOsterloh &: BeckwithI Il995l : [Andrews fc: Williams! l2005l ) , but the detailed distribution of 
-^pair with respect to Op has remained unclear for two reasons: the low resolution of single-dish 
millimeter-wave photometry often included several pairs together, and the multiplicity census of 
the nearest star-forming regions was incomplete. Our survey mitigates these issues, providing an 
opportunity to look at the details of the millimeter luminosity-separation distribution. 

Figure [Ih shows Fpair as a function of ap for the 111 stellar pairs in Taurus (see Table [6]), along 
with the 52 single stars that have millimeter-wave measurements (see Table O black points, gray 
upper limits) . This diagram is a striking confirmation of the original conclusions of Jensen et al. , 
plainly demonstrating that millimeter continuum luminosities scale with the separation between 
stellar pairs. However, as the pair separation decreases, the maximum millimeter luminosities ap- 
pear to decline in discrete jumps (rather than continuously) at two relatively well-defined locations: 
Up sa 300 and 30 AU (p w 2 and 0.2", respectively). These features facilitate a natural breakdown 
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of the Taurus pairs into distinct sub-populations. The distribution of -Fpair for wide pairs, with 
projected separations greater than 300 AU (blue), is similar to the distribution for single stars. At 
medium separations - ap = 30-300 AU (red) - we find a notable absence of bright pairs and a de- 
creased detection rate. At yet smaller separations, Op < 30 AU (green), only a few pairs exhibit very 
weak millimeter-wave emission. A small group of dramatic outliers sparsely populate the otherwise 
empty region of bright pairs with small separa tions (purple): th ese pairs are known or suspected 



to harbor massive ciicuuibinary rings (e.g., see lPietu et al.ll201ll . regarding GG Tau Aab). 



The cumulative distributions of -Fpair foi' each separation-based sub-population can better 
quantify this apparent trend. The 'Pa(>Fpah) distributions shown together in Figure O) were 
constructed using the Kaplan-Meier product-limit estimator to properly accoun t for the substantial 



numb er of pairs in each sub-population that exhibit no millimeter emission (see iFeigelson &: Nelson 



19851 ). The distributions are compared directly with the standard two-sample tests used in survival 
analysis in Table [9l These tests confirm the qualitative examination of Figure [Gh: (1) wide pairs 
and single stars have statistically indistinguishable millimeter luminosity distributions; (2) medium 
pairs have significantly lower luminosities; and (3) small pairs have yet less millimeter emission. 
The 'Pa(>-Fpair) in Figure [6f) have similar functional forms, albeit shifted in luminosity. A simple 
scaling indicates that -Fpair decreases by a factor of ~5 from wide to medium separations, and then 
another factor of 5 from medium to small separations. These trends are not an artifact of including 
non-detections in the analysis: similar conclusions are drawn by comparing only the detected pairs 
in the same sub-populations using the two-sided Kolmogorov-Smirnov (K-S) test (see Table [9]). 

Models of star-disk interactions also postulate an association between the amount of disk trun- 
cation and the component masses of the stellar pair. Massive companions impart larger dynamical 
perturbations to individual disks, producing more tidal stripping, and leading to lower disk masses 
and therefore less millimeter-wave emission. In that case, Fpair should be anti-correlated with q: 
higher mass ratio pairs should have fainter disk emission. No such trend is obvious in Figure [Tji. 
If we separate the full population into high and low mass ratio pairs at some critical Qc, we find 
the largest difference between those sub-categories for qc = 0.5. Figure O) compares the cumulative 
distributions of the millimeter flux densities for high (q > 0.5; blue) and low (q < 0.5; red) mass 
ratio pairs, again with Vq(> -Fpair) constructed using the Kaplan-Meier estimator to include the 
pairs that do not have detectable millimeter emission. The same two-sample tests employed above 
indicate a weak relationship between Fpair and q (see Tabled]), such that pairs with low stellar mass 
ratios have slightly less millimeter emission - the opposite of expectations from tidal interaction 
models. However, the evidence for any increase in the millimeter emission with higher stellar mass 
ratios is contained entirely in the relative detection ratios: a larger fraction of pairs with low q are 
not detected at millimeter wavelengths (see Figure [TJi). If only the detected pairs are compared, 
the millimeter luminosity is found to be independent of q (see the K-S test results in Table [9]). 
Moreover, this trend is present (and in fact enhanced) only for stellar pairs with wide separations: 
no clear relationship between -Fpair and q exists for pairs with ap < 300 AU (see Table [9]) . 

The absence of a firm connection between the millimeter luminosity and stellar mass ratio for 
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Fig. 7. — (a) A comparison of millimeter flux densities from potentially interacting pairs as a 
function of stellar mass ratio. For reference, the pair population is distinguished into high {q > 0.5; 
blue) and low {q < 0.5; red) mass ratio groups. (6) The cumulative distributions of millimeter flux 
densities for high and low mass ratio pairs, constructed with the Kaplan-Meier estimator to include 
the available upper limits. Millimeter luminosity depends only weakly on the stellar mass ratio. 



smaller-separation pairs is consistent with the weak ^-dependence predicted by tidal interaction 
models. Nevertheless, the correlation between -Fpair and q for widely-separated pairs is compelling. 
Since star-disk interactions in these cases should be inherently less destructive, this trend might 
indicate a relationship between the millimeter luminosity and the absolute (and not relative) stellar 
masses in the pair. Indeed, Figure [8]i demonstrates a marginal association between the millimeter 
luminosity and total stellar mass in Taurus pairs. Taking only the detections, we find a ~3 a corre- 
lation between Fpair and Mpair (this improves slightly to 3.7 a if the labeled outliers are excluded). 
Figure [EI) demonstrates that this relationship is not restricted to stellar pairs, but apparently also 
applies to individual stars, both isolated (single) cases and members of multiple systems. In the 
latter case, we have combined the component-resolved F^ estimates from Tableland the single star 
measurements from the literature (Table [7]): the result is again a ~3(T positive correlation (4.5 o" if 
SU Aur is excluded). It is worth noting that the M^, distributions for isolated stars and individual 
stars in multiple systems are statistically indistinguishable. However, it is premature to draw any 
firm conclusions from these weak trends for two reasons: (1) the dispersion is large relative to the 
range of the trend, and (2) the stellar masses used here have large systematic uncertainties. 

The demographic properties of stellar pairs discussed above are certainly informative, but 
they also naturally hide some characteristics of individual components that are available from 
this resolved survey. The two panels in Figure [9] are intended to compare the millimeter- wave 
luminosities from individual disk components within each stellar pair, as a function of their stellar 
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Fig. 8. — (a) Millimeter flux density as a function of the total mass of a stellar pair, where 
-^pair = Mp + Ms = Mp(l + g). (&) Millimeter flux densities as a function of stellar mass for the 
individual, detected stellar components in Taurus, for both single stars (black) and those in multiple 
systems (blue). Note that upper limits for non-detections are not shown in the right-hand panel. 



host masses. Figure [9l3 directly compares the resolved 880 ;um flux densities for each component 
in each pair. As mentioned in §5.1, the emission from the primary is usually more luminous than 
from the secondary. The exceptions above the dashed line are comprised of the widely-separated 
tertiary companions of close pairs and the UZ Tau Wab and FS Tau/Haro 6-5B pairs. The relative 
dominance of the primary disk emission is unaffected by the projected separation to a companion, 
and is typically more than would be expected if the amount of millimeter- wave disk emission scales 
linearly with the stellar host mass. This latter feature is shown more directly in Figure ^, which 
compares the cumulative distributions of the ratio of the millimeter flux density to the stellar mass 
(akin to a disk:star mass ratio) for singles, primaries, and secondaries. In some studies of disks 
around low-mass stars and brown dwarfs, this F^/M^ ratio is found to be roughly constant; the 
weaker emission (or lower detection ra te) is a manifestation of inherently lower host masses (e.g.. 



Scholz et al. 



2006 



Schaefer et al.ll2009l ). That implicit sensitivity threshold is not the case for the 



secondaries in our sample: the Fa/M^, ratio is systematically lower for secondaries compared to 
primaries and singles, with a probability that it is drawn from the same parent distributions as 
the primaries or singles of <0.008 (two-sample test results are compiled in Table [9]). The results 
suggest that the millimeter-wave disk emission from the secondaries is inherently less luminous 
than around the primaries (or isolated stars), regardless of the stellar host mass. 
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Fig. 9. — (a) The 880 /xm flux densities for the primary and secondary components of each stellar 
pair. (6) The cumulative distributions of the ratio of millimeter-wave disk flux densities to their 
stellar host masses (a proxy for the disk:star mass ratio) for singles {black), primaries {blue), and 
secondaries {red), each constructed from the Kaplan- Meier product-limit estimator to include upper 
limits on F^. Two-sample tests (see Table [9]) demonstrate that these ratios are systematically lower 
for secondaries, implying that their disks are inherently less luminous regardless of the host masses. 



5.3. Disk Sizes and Tidal Truncation 



A more direct test of dynamical predictions for star-disk interactions lies with our measure- 
ments of individual disk sizes. Theoretical models provide a way to estimate the truncated equilib- 
rium tidal radii {Rt) for the disks around each component of a stellar pair given a few key orbital 
parameters, ja, e, q|, and some characterization of the viscous properties of the disk material 
(jArtymowicz &: Lubowlll994l ). Although they remain uncertain, there are reasonable ways to es- 
timate stellar mass ratios {q) from optical/infrared measurements and pre-main sequence stellar 
evolution models (see the references in Table ED . However, the stellar pairs in this sample have wide 
enough projected separations that they are expected to have prohibitively long orbital periods and 
exhibit little apparent motion on reasonable time baselines for observations: therefore, we generally 
do not have any direct knowledge about their true orbital separations (a) or eccentricities (e). 

Nevertheless, we can construct a probabilistic model of Rt using only the projected physical 
separation ofany pair, g^ (based on the projected angular separation, p, and assumed distance, d). 
Following [Torrea (119991 ). the ratio of the semimajor axis to the projected physical separation is 



T= — 



1 



1 — sin^ (cj + u) sin^ i, 



(1) 



1 + e cos i^ 
where e is the eccentricity, u is the true anomaly, uj is the longitude of periastron, and i is the orbital 
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inclination relative to the observer (note that the ratio T is exact for binaries; we expect only modest 
deviations from it for well-separated hierarchical pairs). Since we do not know {w, i, e, i^} for any 
individual stellar pair, we have to construct a probability distribution for the true-to-projected 
separation ratios, V{J-\ using a Monte Carlo approach. To accomplish that, we assume that 
stellar pairs are not observed at any preferential orbital location and adopt uniform distributions 
for the orbital phase (or mean anomaly) and longitude of periastron (w). The assumption of random 
viewing geometries suggests that the distribution of orbital inclinations {%) has a sin i dependence. 
However, inferring an appropriate functional form for the eccentricity (e) distribution (and by 
extension the distribution of true anomalies, z^) is more challenging. 

There is little empirical information available to constrain the eccentricity distribution for 
the pre-main sequence binary population. Pairs with short or bital period s have low eccentricities 



1977 



Zahn &: Bouchet 



198S 



(e < 0.1), due to the rapid tidal circularization of their orbits (jZahn 

Melo et al.ll200ll ). At longer period s, the eccentricity distribution appears relatively uniform in 

the range e ~ 0.1-0.9 (JMathieulll994l ). It is not clear if the apparent dearth of young stellar pairs 



with extreme (circular or parabolic) eccentricities is real or an artifact of low-number statistics 
and se lectio n effects. Based o n their samples of main sequence field binaries, iDuquennov &: Mavor 
(jl99ll ) ai id iTokovininI (119981) s uggested an increasing eccentricity distribution, where V{e) oc 2e 



(see also 



Ambartsumian 



19371 ). However, the recent comprehensive survey of such systems by 



Raghavan et al.l (|2010l ) instead suggests that a flat eccentricity distribution is prefer able for orbital 
perio ds of ~10-10^ days. Similar results are also noted for very low mass binaries (JDupuv &: Liu 
201 ll ). Based on these more recent studies, we assume that the eccentricity distribution is uniform. 
For reference. Figure [TOh illustrates the influence of different forms for V{e) on the shape of V{J-). 

Having established the infrastructure to derive a probabilistic model of the true orbital sep- 
aration (a) given the observed projected separation (flp), we can build on that to determine the 
distribution of tidal truncation radii, V{Rt), for a given stellar pair. Because of its relative sim- 
plicity, we adopt the sem i-analytic approximati ons for truncated disk sizes based on the analysis of 
stable invariant loops by lPichardo et al.l (120051 ). Using their formulation, the tidal radius is 



Rt ^ 0.337 



(1 



g)1.20 ^2/3 ^0.07 



0.6c^2/3+ln(l + (^l/3) 



J-" a 



P 1 



(2) 



where /x = q/i^ + q) is the mass fraction of the stellar pair and 99 is the mass ratio of the host 
star for which Rt is being calculated relative to its companion. For example, the truncation radius 
of the disk around the primary star is calculated with ip = Mp/Mg = l/g, whereas the Rt for the 
disk around the secondary is determined by setting ip = Mg/Mp = q. For any pair in a multiple 
system, there is a direct measurement of Op (or rather p, and an assumed distance) and an estimate 
of q (see Table [6]). Fixing these quantities, we constructed the probability distribution V{Rt) that 
a component of the pair hosts a disk with a tidally truncated radius Rt using a Monte Carlo 
simulation with ~10^ realizations of Equation (2), assuming the priors for the distributions of the 
orbital elements {w, i, e, z^} (and therefore J-) discussed above. 
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Fig. 10. — (a) Probability distributions for the ratio of true to projected separations, V{J^) where 
J^ = a/ttp, using three different underlying eccentricity distributions. (6) An example probability 
distributi on for the disk radii (R t) expected in the HK Tau binary based on the tidal interaction 
models of iPichardo et al.l (|2005l ) for the same three assumed eccentricity distributions. The blue 
and red vertical bars represent the best-fit disk radii (and 1 a uncertainties) for HK Tau A and B, 
respectively, measured directly from the SMA visibilities in §4. The yellow filled area marks the 
region containing 68% of the probability for the V{Rt) with the favored eccentricity distribution 
(black curve; see text), used to determine the error bars in Figure [TTl 



As an example. Figure [TOb shows V{Rt) for the HK Tau binary, where Op = 340 AU and (7 ~ 1, 
for three representative assumptions about the underlying eccentricity distribution (note that the 
same Rt distribution applies to both components for this equal- mass stellar pair). The best-fit 
estimates of the disk radii from our modeling of the SMA data (see §4) are marked as red (HK 
Tau A) and blue (HK Tau B) vertical bars. Tidal interaction models predict that the disk sizes 
have a rather steep dependence on the orbital eccentricity (see Equation 2), which means that the 
assumption of an underlying V{e) that permits or favors high eccentricities will lead to the general 
prediction of very small disk sizes (gray or green curves in FigureflOb) and, therefore, low millimeter- 
wave luminosities. While such eccenticity distributions may be relevant for the general population 
of multiple systems, the luminosity-based selection criterion used to build our component-resolved 
SMA sample creates a strong bias that would exclude high-e pairs. With that bias in mind, we 
favor the use of a truncated eccentricity distribution to make comparisons between the measured 
and predicted disk radii; we assume V{e) is uniform for e G [0.0,0.7] (black curves in Figure [TOj) . 

Figure [TT] makes a direct comparison between the measurements of dust disk sizes {Rd) that 
were determined in §4 (s ee Tables HI and [ p and the truncation radii {Rt) predicted by our prob- 
abilistic treatment of the IPichardo et al.l (J2005l ) models. The location of the points along the Rt 
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Fig. 11. — The measured disk radii (see Tables [H and [5l) compa red with the expected disk radii 
based on a probabilistic treatment of the iPichardo et al.l (|2005l ) tidal interaction models, for an 
assumed underlying uniform eccentricity distribution truncated at e = 0.7 (see text). 



axis (abscissae) correspond to the peaks of their V{Rt) distributions (for the assumed V{e) de- 
scribed above), and their asymmetric error bars encapsulate the central 68% of those probabihty 
distributions (see the shaded yeUow region of Figure [TOb for an example). Of the 14 disks with 
available Rd measurements, eight have sizes that are in good agreement with predictions from the 
tidal interaction models, two are smaller than expected (around GG Tau Aab and DK Tau A), and 



the re maining four are considerably larger. These decidedly mixed results for the iPichardo et al 



(J2005l ) model predictions could be at least partially ameliorated by introducing a term that in- 
corporates viscosity into the interaction calculations, which would t end to increase Rt due to the 
viscous spreading of disk material (see lArtvmowicz fc Lubowl Il994l ) . Quantitatively comparable 
shifts in the measured Rd values could be accomodated by permitting a range of emission gradients 
(x) in model fits to the observations. Reconciling the measured disk sizes with the predictions from 
tidal interaction models would require that the disk viscosity increases and/or the millimeter- wave 
emission gradient (x) decreases as a function of pair separation (a). 



6. Discussion 



We have carried out a luminosity and separation limited survey of the millimeter-wave dust 
continuum emission from the disks around individual components of 23 young multiple star systems 
in the ~l-2 Myr-old Taurus- Auriga star formation region. With a simple morphological model, we 
fitted the continuum visibilities observed with the SMA interferometer to determine the luminosity 
and size of each individual disk in these systems. These component-resolved measurements were 
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then coupled with a comprehensive database of milhnieter-wave luminosities for all of the multi- 
ple systems in Taurus (with spectral types F0-M4) to estimate millimeter detection frequencies, 
evaluate the dependence of continuum emission levels on the separations and masses of companion 
stars, and make direct comparisons with predictions from tidal interaction models. 

We find that roughly one third (28-37%) of the individual stars in multiple systems harbor 
disks with dust masses that are large enough to emit detectable millimeter continuum radiation 
(see §5.1). This low incidence rate is approximately half that for isolated single stars (62%), and 
does not depend on the number of stellar companions in the sys tem. Similar disk frequencies 



have been inferred from accretion and infrared excess signatures (|Cieza et al.l l2009l : iKraus et al 



2012I ). although without component-resolved diagnostics. Given these disk frequencies, it is clear 
that some external processes driven by the presence of a companion act to hasten the dispersal 
of circumstellar material in multiple systems at a level comparable to any internal disk evolution 
mechanisms (e.g., photoevaporation, grain growth, planet formation, etc.) on ~l-2Myr timescales. 

Some basic demographic properties of stellar pairs in Taurus can provide new insights into the 



mechanics of those external evolution processes. Building on the initial work bv lJensen et al.l (|1994| . 



1996bl ). we have shown that the millimeter- wave luminosity from a pair of stars depends strongly 
on their (projected) separation (see §5.2). We identified substantial changes in the millimeter lumi- 
nosity distributions of pair populations at discrete separations of 30 and 300 AU. Widely-separated 
pairs (cp > 300 AU) have emission levels similar to single stars. Pairs with medium separations 
{ttp = 30-300 AU) are typically 5x fainter, with a lower overall detection fraction. Millimeter emis- 
sion is rarely detected around pairs with small separations (op < 30 AU), representing at least 
another factor of ~5 reduction in luminosity. We demonstrated that there is a weak tendency 
for pairs with comparable stellar masses (higher q) to have brighter millimeter emission, with the 
effect being considerably stronger for wider pairs. We suggested this is related to a marginal cor- 
relation between stellar mass and millimeter luminosity, although verifying that tentative trend is 
a challenge due to the systematic uncertainties involved in estimating stellar masses. 

The relationship between millimeter luminosity and pair separation (Fig. [6]) suggests that the 
external process relevant for the evolution of circumstellar material in multiple systems may be tie d 
to tidal stripping from the outer regions of their constituent disks (e.g.. lArtvmowicz Sz Lubowlll994l ). 



However, we found mixed results in a comparison of individual disk sizes that were measured (from 
the SMA data; §4) and predicted from tidal interaction models (see §5.3). Using a probabilistic 
treatment of orbital parameters, we showed that about half of the resolved dis k s in o ur sample have 



sizes that are consistent with the truncation predictions of the lPichardo et al.l (|2005l ) models. Most 
of the remainder have sizes that are substantially larger than expected, given the smaller projected 
separations to their companions (Fig. [TTI) . Analogous discrepan cies have been noted regarding the 



obser ved and predicted inner edges of circumbinary rings (e.g.. lBeust &: Dutreyll2005l : lNagel et al 



20101 ) . These results hint that another external process also shapes the circumstellar environments 
of young multiples. Further support for this additional evolutionary process is found in the observed 
distributions of circumstellar dust around individual stellar components. Tidal interaction models 
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predict that disks should survive around both components of a pair, with their relative sizes roughly 
set by q (for a given a). However, Fig. [9] (see §5.2) demonstrates that this is often not the case. 
Roughly half of the multiple systems with detectable millimeter emission harbor only a single disk, 
usually around the primary component (or the wide tertiary in some hierarchical triples; e.g., HV 
Tau or Haro 6-37). Moreover, aside from a few exceptional counterexamples (see Fig. [Gji), we find 
little millimeter-wave evidence for the circumbinary disks that should be common around pairs 
with small separations if tidal interactions were the sole external evolution mechanism. 

To be fair, tidal interactions alone may still be able to explain many of these observed proper- 
ties. The processes of stripping and truncation in pairs where the orbital and disk planes are mis- 
aligned has not yet been explored in detail (although see, e.g.. lAkeson et al.ll2007l : IVerrier &: Evans 
20081 ). but might be substantially enhanced for some configurations. With limited orbital infor- 
mation, the prevalence of such misalignment is not known: bu t, there is some indication from 
polarization measurements that it is common (j Jensen et al.l l2004ll and a nu i nber of specific exarn - 



ples have been identified (e.g., HK Tau or HD 98800; iDuchene et al.ll200,i lAndrews et al.ll2010l V 
Moreover, it is important to keep in mind that the emission we have measured is only capable of 
probing the trace population of dust particles in these disks, and not their dominant mass reser- 
voirs of molecular gas. The total dis k masses that would be in f erred f rom this emission could 



be substantiall y under-estimated (see lAndrews &: Williams 



disk sizes (e.g.. lHughes et al.ll2008l : IPanic et al. 
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2005, 



Andrews et al. 



2007a '). as could the apparent 



2OI2I ). Future complementary 



observations of emission line probes would help better address such uncertainties. 

Alternatively, the disk properties we observe may be set at very early stages by the processes 
that regulate accretion during multiple star formation. Numerical simulations indicate that the ratio 
of the specific angular momenta of the infalling material (jgas) to the stellar pair (j^,) is fundamental 
in determining how the gas an d dust accreted from a proto-system envelop e is distributed among 



individual stellar components (JBonnell fc Batdll994l : iBate fc Bonnelll 119971 ). If jgas < j*, rnost of 



the infalling material will form a disk around the stellar primar y. Conversely, ii jpas > j*, then 
a circum-secondary or circumbinary disk will dominate (see also lOchi et al.l l2005l ) . In the many 
cases with only circum-primary disk detections described above, the data are consistent with the 
former scenario. Similar millimeter-wave observations suggest that primary stars preferentially 
harb or more circumstellar material at the even earlier Class O/I stages of protostellar evolution 



(e.g.. lLaunhardtll2004l : iPatience et al.ll2008l . and references therein). 



To summarize, there is a body of observational evidence suggesting that (at least) two funda- 
mental processes related to the presence of stellar companions play significant roles in the evolution 
of circumstellar material in young multiple systems. The first is associated with the multiple star 
formation process itself, where the fraction of angular momentum associated with infalling mate- 
rial relative to that contained in the orbital motion of the stellar pair determines how circumstellar 
material is apportioned to each component. This is likely responsible for the pairs we observe with 
very large primary-to-secondary millimeter luminosity ratios. The second is a tidal interaction 
process that strips material from any individual disks that survive the formation process. This is 
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thought to be the origin of the minimeter luminosity - pair separation relationship. 

At the typical ~l-2Myr age of Taurus- Auriga multiples, both of these evolution mechanisms 
have already made their mark on disk properties. The long-term fate of the circumstellar material 
in these multiple systems now rests with the same internal mechanisms that govern the subsequent 
evolution and dissipation of the disks around their isolated (single) counterparts. The formation 
and evolution of planets from this disk material are mechanisms of particular interest in these 
systems. The relatively straightforward, albeit uncertain, relationship between millimeter-wave lu- 
minosities and total disk masses enables at least a preliminary assessment of the likelihood of planet 
formation around ind i vidual comp onents in multiple systems. Using the conversion advocated by 
Andrews &: Williamsl (|2005l . l2007al ) - which assumes optically thin, isothermal dust emission with 



a mean temperature of 20 K and (gas-|-dust) opacity of 0.034 cm^ g~^ at 880 ;um - Fd ^ 15mJy 
corresponds to a disk mass of ~lMjup, and F^ ~ 150 mJy represents the standard estimate of the 
minimum mass of the solar nebula, ~0.01Mq. Taking these conversions at face value, we would 
conclude that giant planets are unlikely to form around stars in close pairs (Fig. [6]) or around the 
secondary components of most pairs with wider separations (Fig. [9]). However, we would likewise 
infer that the primary components in wider pairs, the wide tertiaries in hierarchical triples, and 
perhaps the population of spectroscopic binaries should be just as likely to host giant planets as 
single stars. 

Direct exoplanet searches of main sequence multiple systems in the field confirm these gen- 
eral expectations from the disk survey: planet formation is not severely inhibited by the pres- 



ence of a stellar companion (e.g., iPatience et al. 



2002 



Raghavan et al.l 12006 : IPuchend l20ld ) , and 



giant plane ts are preferentially found around primaries, the isolated components of hierarchi 



cal triples (IDesidera &: Barbieri 



Desidera et al. 



2011 



2007 



Bonavita &: Desidera 



2007 



Mugrauer &: Neuhauser 



I, and perhaps even spectroscopic binaries (e.g.. lDovle et al.ll201ll ). That cor- 
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roboration of results is a promising sign for planet formation studies in multiple systems, but the 
characterization of the disks around individual stellar components in these systems is still in its 
early stages. In the near future, we expect that new facilities like the Atacama Large Millime- 
ter Array (ALMA) will shift the focus to study these individual disks in detail to directly compare 
their density structures and particle growth signatures with the disks around isolated (single) stars. 
Ultimately, such observations will help determine the impact of a nearby stellar neighbor on the 
planet formation process. 



7. Summary 



We have presented the results of an SMA imaging survey of the millimeter-wave (880 /xm or 
1.3 mm) thermal continuum emission from circumstellar dust in 23 young multiple systems in the 
Taurus- Auriga star-forming region. This census was designed to target relatively bright (>20mJy 
at 880 ;um) and well-separated (p > 0'.'3) systems with primary spectral types between FO and M4. 
We employed simple morphological models of the SMA visibilities to measure component-resolved 
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millimeter luminosities and disk sizes whenever possible. Those results were considered together 
with a comprehensive literature compilation of the millimeter luminosities from all Taurus multiples 
(in this spectral type range) to better analyze how the presence of a stellar compaion affects basic 
disk properties. Our primary conclusions are the following: 

1. The millimeter detection frequency for individual stars in multiple systems is approximately 
1/3 (28-40%), about half that for single stars (62%), and independent of the number of 
companions. These relative incidence rates suggest that the presence of a stellar companion 
plays a substantial external role in the early development and evolution of circumstellar 
material, at a level comparable to the standard internal disk evolution mechanisms that can 
operate in isolation (e.g., photoevaporation, particle growth, planet formation, etc.). 

2. The millimeter luminosity from a pair of stars depends strongly on their projected separation 
(op), such that closer pairs are substantially fainter. We find natural breaks in the luminosity- 
separation plane at ap ^ 30 and 300 AU. The luminosity distribution of wide pairs (op > 
300 AU) is indistinguishable from that of single stars. Pairs with medium separations (cp = 
30-300 AU) are 5x fainter, and the very few close pairs (op < 30 AU) that we detect are 5x 
fainter yet - although a few bright circum6mary disks represent notable exceptions. 

3. There is no clear relationship between the millimeter luminosity from a pair and its stellar 
mass ratio (g) in general, but wide pairs with higher q tend to be brighter. We show that this 
latter behavior is produced by a marginal correlation between millimeter luminosities and 
stellar masses (both summed among pairs and for individual stars). However, the significance 
of this trend is questionable: the scatter is large, and there are substantial systematic un- 
certainties in estimating individual stellar masses that are not considered here. In nearly all 
cases, the primary component of a binary pair or the wide tertiary of a hierarchical triple har- 
bors the disk material that dominates the millimeter luminosity of the system (higher-order 
systems show a range of behavior, depending on the hierarchical nature of their pairings). 

4. We find mixed results from a direct comparison of the disk sizes measured from the data 
and predicted from tidal interaction models, based on a probabilistic treatment of the orbital 
parameters for each syst em. Of the 15 resolve d disks in our sample, the radii expected from 



the models described bv lPichardo et al.l (|2005l ) are commensurate with the observed radii for 



eight disks; two others are found to be too small, and the remaining five are notably larger 
than would be expected given the relatively small (projected) distances to their companions. 

5. These millimeter- wave observations suggest that at least two external mechanisms contribute 
to the evolution of circumstellar material in young multiple star systems. Star-disk tidal 
interactions strip material from the outer regions of individual disks, a process responsible 
for the strong dependence of the millimeter luminosity on the projected separation to a 
companion. We are lead to infer that accretion during the multiple star formation process 
itself also plays a substantial role in apportioning disk material, in many cases setting the 
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initial disk masses such that primary stars harbor substantially more circumstellar material 
than their companions. The long-term prospects for (giant) planet formation in multiple 
systems are poor for stars in close pairings and around secondaries, but should be comparable 
to those around single stars for primaries, wide tertiaries in hierarchical triples, and perhaps 
around both components of very close spectroscopic binaries. 



We are grateful to Trent Dupuy for his assistance with the probabilistic treatment of projected 
orbits, to Joanna Brown for kindly providing some supplementary observing time on the FQ Tau 
binary, and to an anonymous referee for helpful suggestions on clarifying the draft manuscript. The 
Submillimeter Array (SMA) is a joint project between the Smithsonian Astrophysical Observatory 
and the Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 
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Table 1. Selected Sample of Taurus Multiples 



System 


iV. 


i Vpair 


MHO 1/MHO 2 


3 


2 


DD Tau 


2 


1 


FQ Tau 


2 


1 


TTau 


3 


2 


FS Tau/Haro 6-5B 


3 


2 


FVTau 


4 


3 


DI Tau/DH Tau 


4 


3 


UXTau 


4 


4 


FXTau 


2 


1 


DK Tau 


2 


1 


XZ Tau/HL Tau 


3 


2 


HKTau 


2 


1 


V710 Tau 


3 


2 


GGTau 


4 


3 


FZ Tau/ FY Tau 


2 


1 


UZTau 


4 


3 


JH 112 


3 


2 


V807 Tau/GH Tau 


5 


4 


GI Tau/GK Tau 


2 


1 


HNTau 


2 


1 


IT Tau 


2 


1 


HP Tau 


4 


4 


HVTau 


3 


2 


Haro 6-37 


3 


2 


UY Aur 


2 


1 


CIDA-9 


2 


1 


RW Aur 


2 


1 



Note. — Col. (1): System name. Italicized systems 
have not been re-observed with the SMA (see §2). 
Col (2): Number of stars in system. Col. (3): Number 
of "pairs" in system (see §2). 
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Table 2. SMA Observing Log 



UT Date Targets Config. i/(LO) 

[GHz] 

2005 Dec 17 DH/DI Tau 

2009 Oct 24 DK Tau 

2010 Feb 11 DK Tau 
2010 Feb 18 UZ Tau 
2010 Feb 19 HK Tau 
2010 Mar 2 UZ Tau 
2010 Nov 10 IT Tau, FX Tau, HK Tau, 

V710 Tau, HN Tau, CIDA-9 
2010 Nov 11 FQ Tau, IT Tau, FV Tau, 

HV Tau, JH 112 
2010 Nov 12 FQ Tau, IT Tau, FV Tau, 

HV Tau, JH 112 
2010 Dec 13 RW Aur, Haro 6-37, DD Tau, 

T Tau, UY Aur 
2010 Dec 25 HP Tau 

2010 Dec 27 GG Tau, GI/GK Tau, MHO 1/2 

2011 Jan 28 FQ Tau 
2011 Feb 3 FV Tau, FX Tau, T Tau 
2011 Feb 17 UZ Tau, UY Aur, DD Tau 

2011 Aug 16 DD Tau 

2012 Jan 24 FV Tau, FX Tau, UY Aur 
2011 Oct 13 J04080782+2807280, MHO 3 C 225.497 

J04141 188+281 1535, V410 X-ray 7 
2011 Oct 27 V710 Tau C, J04414565-f2301580 C 222.709 

J04554757+3028077, CFHT Tau 7 



Note. — Col. (1): UT date of observations. Col. (2): Multiple systems 
observed. Col. (3): SMA antenna configuration (see §3). Col. (4): LO 
frequency setting. Note that the observations below the horizontal line 
focused on multiple systems that are not part of our sample: they were 
observed because there were no previous millimeter-wave measurements of 
these systems available in the literature (see §5.1). 



c 


340.755 


c 


340.767 


V 


343.277 


V 


341.536 


V 


342.289 


V 


340.427 


c 


340.755 


c 


340.755 


c 


340.755 


c 


340.755 


c 


225.497 


c 


225.497 


E 


340.755 


E 


340.755 


E 


340.755 


E 


341.575 


E 


340.755 
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Table 3. Properties of Synthesized Continuum Maps 



System 


A 


RMS noise 


beam size 


beam PA 




[mH 


[mjy beam^ ] 


["] 


n 


(1) 


(2) 


' (3) 


(4) 


(5) 


MHO 1/2 


1335 


2.1 


2.8 X 2.2 


95 


DD Tau 


880 


1.9 


0.9 X 0.7 


79 


FQ Tau 


880 


1.2 


0.9 X 0.7 


77 


T Tau 


880 


3.6 


1.0 X 0.8 


77 


FVTau 


880 


1.1 


1.0 X 0.7 


79 


DI/DH Tau 


880 


3.3 


2.6 X 2.1 


151 


UXTau 


880 


0.9 


0.6 X 0.4 


33 


FXTau 


880 


1.6 


1.1 X 0.8 


82 


DK Tau 


880 


0.7 


0.5 X 0.5 





HK Tau 


880 


1.4 


0.7 X 0.6 


65 


V710 Tau 


880 


3.2 


2.2 X 2.0 


94 


GG Tau 


1335 


2.2 


3.2 X 2.4 


32 


UZTau 


880 


0.7 


0.4 X 0.3 


51 


JH 112 


880 


1.9 


2.0 X 1.8 


57 


GI/GK Tau 


1335 


0.7 


3.6 X 2.6 


68 


HNTau 


880 


2.1 


2.2 X 2.0 


96 


IT Tau 


880 


1.4 


2.1 X 1.9 


78 


HP Tau 


1335 


1.2 


3.9 X 2.9 


102 


HVTau 


880 


2.3 


2.0 X 1.8 


61 


Haro 6-37 


880 


5.0 


2.1 X 2.0 


136 


UY Aur 


880 


1.5 


0.9 X 0.6 


91 


CIDA-9 


880 


2.4 


2.2 X 2.0 


98 


RW Aur 


880 


3.6 


2.1 X 2.0 


159 



Note. — Col. (1): System name. Col. (2): Effective observing 
wavelength. Col. (3): RMS noise level (Icr) in the naturally 
weighted synthesis image (see Figures [2]|4)). Col. (4) & (5): Beam 
dimensions and major axis position angle. 
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Table 4. Resolved Disk Properties 



Name 


Disk Center 


Fi 


Rd 


id 


PAd 




a [J2000] 


5 [J2000] 


[mJy] 


[AU] 


n 


n 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


MHO 1 


04 14 26.28 


+28 06 02.90 


(242 ± 3) 


131 + 28 


32 + 15 


84 + 28 


T TauN 


04 21 59.44 


+19 32 06.31 


304 + 3 


64 + 3 


6 + 21 


90 + 36 


UX Tau A 


04 30 04.00 


+18 13 49.25 


141 ±1 


57+1 


34 + 2 


174 + 4 


FX Tau A 


04 30 29.64 


+24 26 44.77 


24 + 3 


219 + 57 


[45] 


[90] 


DK Tau A 


04 30 44.25 


+26 01 24.52 


57 + 1 


76 + 5 


27 + 9 


115 + 12 


HK Tau A 


04 31 50.57 


+24 24 17.60 


81 ±2 


102 + 8 


57 + 4 


118 + 5 


HK Tau B 


04 31 50.60 


+24 24 15.27 


49 + 2 


164 + 22 


[84] 


[40] 


GG Tau Aab 


04 32 30.36 


+17 31 40.03 


(557 ± 3) 


305 + 5 


[37] 


[7] 


UZ Tau Eab 


04 32 43.07 


+25 52 30.80 


367 + 3 


219 + 3 


50 + 1 


86 + 1 


UZ Tau Wa 


04 32 42.83 


+25 52 31.04 


24 + 2 


63 + 6 


56 + 4 


40 + 6 


UZ Tau Wb 


04 32 42.83 


+25 52 31.41 


41 + 2 


119 + 10 


[45] 


[90] 


HN Tau A 


04 33 39.38 


+17 51 52.12 


34 + 3 


233 + 63 


[45] 


[90] 


HV Tau C 


04 38 35.47 


+26 10 41.24 


98 + 3 


139 + 21 


[80] 


99 + 9 


Haro 6-37 B 


04 46 59.09 


+17 02 39.87 


218 + 6 


219 + 15 


[45] 


[90] 



Note. — Col. (1): Stellar component name. Col. (2) and (3): Adopted disk center 
coordinates. Col. (4): Flux density. Values in parenthesis correspond to measurements at 
a wavelength of 1.3mm: all others are at 880 /u.m. Col. (5): Disk radius. Col. (6): Disk 
inclination angle, where 0° is face-on. Col. (7): Position angle (E of N) of the disk major 
axis projected on the sky. Values in square brackets were fixed in the modeling (see §4). 



Table 5. Unresolved Disk Properties 



Name 


Disk Center 


Fd 


Name 


Disk Center 


Fd 




a [J2000] 


5 [J2000] 


[mJy] 




a [J2000] 


5 [J2000] 


[mJy] 


MHO 2 AB 


04 14 26.47 


+28 05 59.95 


(138 + 1) 


GG Tau Bb 


04 32 30.51 


+17 31 28.68 


(<7) 


DD Tau A 


04 18 31.13 


+28 16 28.82 


24 + 3 


FZ Tau 


04 32 31.76 


+24 20 03.10 


29 + 7 


DD Tau B 


04 18 31.13 


+28 16 28.26 


9 + 3 


FY Tau 


04 32 30.58 


+24 19 57.29 


< 27 


FQ Tau A 


04 19 12.81 


+28 29 32.75 


8 + 2 


JH 112 Aa 


04 32 49.12 


+22 53 02.52 


10 + 3 


FQ Tau B 


04 19 12.79 


+28 29 33.44 


< 5 


JH 112 Ab 


04 32 49.22 


+22 53 03.04 


5 + 3 


T Tau Sab 


04 21 59.44 


+19 32 05.61 


< 10 


JH 112 B 


04 32 49.27 


+22 53 07.89 


<6 


FS Tau AB 


04 22 02.18 


+26 57 30.49 


49 + 6 


V807 Tau AB 


04 33 06.64 


+24 09 54.90 


20 + 3 


Haro 6-5B 


04 22 00.70 


+26 57 32.50 


(134 ± 6) 


GH Tau AB 


04 33 06.22 


+24 09 33.99 


15 + 3 


FV Tau A 


04 26 53.53 


+26 06 54.10 


15 + 1 


GK Tau 


04 33 34.57 


+24 21 05.64 


(3 + 1) 


FV Tau B 


04 26 53.47 


+26 06 54.12 


11 + 1 


GI Tau 


04 33 34.06 


+24 21 16.83 


(12 + 1) 


FV Tau/c A 


04 26 54.41 


+26 06 50.77 


< 5 


HN Tau B 


04 33 39.24 


+17 51 49.71 


<6 


FV Tau/c B 


04 26 54.36 


+26 06 51.03 


< 5 


IT Tau A 


04 33 54.71 


+26 13 27.36 


16 + 2 


DI Tau AB 


04 29 42.46 


+26 32 49.38 


< 10 


IT Tau B 


04 33 54.59 


+26 13 25.21 


11 + 3 


DH Tau A 


04 29 41.64 


+26 32 58.34 


47 + 4 


HP Tau/G3 AB 


04 35 53.50 


+22 54 09.84 


(<4) 


DH Tau B 


04 29 41.74 


+26 32 56.56 


< 10 


HP Tau/G2 


04 35 54.16 


+22 54 13.48 


(<4) 


UX Tau Bab 


04 30 03.59 


+18 13 49.18 


<3 


HP Tau 


04 35 52.77 


+22 54 23.17 


(48 ± 2) 


UX Tau C 


04 30 04.00 


+18 13 43.39 


<3 


HV Tau AB 


04 38 35.29 


+26 10 38.16 


< 7 


FX Tau B 


04 30 29.58 


+24 26 45.06 


< 5 


Haro 6-37 Aab 


04 46 58.98 


+17 02 37.90 


< 15 


DK Tau B 


04 30 44.39 


+26 01 23.45 


<3 


UY Aur A 


04 51 47.39 


+30 47 13.31 


39 + 3 


XZ Tau AB 


04 31 40.07 


+18 13 57.18 


(14 + 3) 


UY Aur B 


04 51 47.34 


+30 47 12.71 


10 + 3 


HL Tau 


04 31 38.44 


+18 13 57.65 


2360 ± 90 


CIDA-9 A 


05 05 22.82 


+25 31 30.65 


77 + 3 


V710 Tau A 


04 31 57.80 


+18 21 37.78 


125 + 3 


CIDA-9 B 


05 05 22.97 


+25 31 31.84 


<7 


V710 Tau B 


04 31 57.81 


+18 21 34.61 


< 10 


RW Aur A 


05 06 49.57 


+30 24 04.90 


62 + 6 


V710 Tau C 


04 31 59.68 


+18 21 30.50 


< 22 


RW Aur B 


05 06 49.47 


+30 24 04.53 


9 + 5 


GG Tau Ba 


04 32 30.30 


+17 31 29.69 


(<7) 











00 

o 



Note. — Col. (1): Stellar component name. Col. (2) and (3): Disk center coordinates. Col. (4): 8 80 ^m flux density. Values i n 
parenthesis correspond to 1.3mm. The sources in italics have Fd values from single-dish photometry ( Andrews fc Williamsll2005h . 
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Table 6. Properties of Stellar Pairs in Taurus 



System 




Pair 


P 


q 


Mp 


-impair 


Refs 




Primary 


Secondary 


["] 




[M0I 


[mJy] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


J04080782+2807280 


J04080782 A 


J04080782 B 


0.044 


0.13 


0.30 


(<8) 


1, 1, A 


V773 Tau 


V773 Tau Aa 


V773 Tau Ab 


0.0028 


0.86 


1.54 


<9 


2, 2, B^ 




V773 Tau Aab 


V773 Tau B 


0.14 


0.29 


2.87 


<9 


3, 3*= ,B^ 




V773 Tau Aab 


V773 Tau C 


0.21 




2.87 


<9 


3, 3*= .B'" 




V773 Tau B 


V773 Tau C 


0.20 




1.33 


<9 


3^^ ,3'" ,B=- 




V773 Tau ABC 


J0414118+2811535 


23.38 


-0.03 


>4.20 


(<9) 


4, I'' ,A 


Anon 1 


Anon 1 A 


Anon 1 B 


0.015 


0.68 


0.64 


<8 


1, 1, B 


MHO 1/2 


MHO 1 


MHO 2 AB 


3.93 


1.00 


0.44 


(380 ± 3) 


4, 1, A 




MHO 2 A 


MHO 2B 


0.050 


0.34 


0.33 


(138 ±1) 


1, 1, A 


MHO 3 


MHO 3 A 


MH0 3B 


0.031 


0.81 


0.72 


«12) 


1, 1, A 


LkCaS 


LkCa 3 Aa 


LkCa 3 Ab 


0.0004 


0.58 


0.36 


<9 


5" .S'^'^B 




LkCa 3 Aab 


LkCa 3 B 


0.48 


0.81 


0.57 


<9 


1, 1, B 


FO Tau 


FO Tau A 


FO Tau B 


0.15 


1.00 


0.33 


13 ±3 


1, 1, B 


LkCa5 


LkCa 5 A 


LkCa 5 B 


0.048 


0.09 


0.50 


<4 


1, 1, B 


V410 Tau 


V410 Tau A 


V410 Tau B 


0.12 


0.20 


0.94 


<8 


1, 1, B^ 




V410 Tau AB 


V410 Tau C 


0.29 


0.08 


1.13 


<8 


1, 1, B" 


DD Tau 


DD Tau A 


DD Tau B 


0.55 


1.00 


0.33 


33 ±4 


1, 1, A 


CZ Tau 


CZ Tau A 


CZ Tau B 


0.32 


0.50 


0.54 


<9 


1, 1, B 


V410 Tau X-ray 7 


V410 X-ray 7 A 


V410 X-ray 7 B 


0.032 


0.65 


0.60 


(<13) 


1, 1, A 


Hubble 4 


Hubble 4 A 


Hubble 4 B 


0.028 


0.73 


0.72 


<9 


1, 1, B 


FQ Tau 


FQ Tau A 


FQ Tau B 


0.75 


0.83 


0.40 


8±2 


1, 1, A 


LkCa7 


LkCa 7 A 


LkCa 7 B 


1.02 


0.52 


0.64 


<9 


1, 1, B 


TTau 


T Tau Sa 


T Tau Sb 


0.082 


0.22 


2.73 


< 10 


6, 6, A" 




T Tau Sab 


TTauN 


0.68 


0.63 


3.34 


304 ±3 


6"* ,6'' ,A 


FS Tau 


FS Tau A 


FS Tau B 


0.23 


0.52 


0.64 


49 ±6 


1, 1, B 




FS Tau AB 


Haro 6-5B 


19.88 


0.85 


0.97 


(155 ±7) 


4, 1, B/C 


LkCa 21 


LkCa 21 A 


LkCa 21 B 


0.044 


0.77 


0.41 


< 10 


1, 1, B 


Jl-4872 


Jl-4872 Aa 


Jl-4872 Ab 


0.17 


0.86 


0.64 


<8 


1, 1, B 




Jl-4872 Ba 


Jl-4872 Bb 


0.10 


0.67 


0.59 


<8 


7, 7'' ,B 




Jl-4872 Aab 


Jl-4872 Bab 


3.38 


0.83 


1.19 


<8 


4, 1, B 


FVTau 


FV Tau A 


FV Tau B 


0.70 


0.73 


0.82 


26 ±2 


1, 1, A 




FV Tau/c A 


FV Tau/c B 


0.70 


0.73 


0.45 


<5 


1, 1, A 




FV Tau AB 


FV Tau/c AB 


12.29 


0.55 


1.42 


26 ±2 


4, 1, A 


DF Tau 


DF Tau A 


DF Tau B 


0.073 


0.90 


0.50 


9±2 


1, 1, B 


Jl-507 


Jl-507 A 


Jl-507 B 


0.079 


0.98 


0.27 


<6 


1, 1, B 


FWTau 


FW Tau A 


FW Tau B 


0.15 


1.00 


0.27 


5±1 


1, 1, B 


DI/DH Tau 


DI Tau A 


DI Tau B 


0.12 


0.13 


0.64 


< 10 


1, 1, A 




DH Tau A 


DH Tau B 


2.34 


0.07 


0.57 


47 ±4 


8, 8*^ ,A 




DI Tau AB 


DH Tau AB 


15.23 


0.85 


0.72 


47 ±4 


4, 1, A 
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Table 6 — Continued 



System 


P< 


iir 


P 


<? 


Mp 


J^pair 


Refs 




Primary 


Secondary 


n 




[Mq] 


[mJy] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


UXTau 


UX Tau A 


UX Tau C 


2.69 


0.15 


1.20 


141 ±1 


1, 1, A 




UX Tau A 


UX Tau Bab 


5.86 


0.86 


1.20 


141 ±1 


1, 1 


A 




UX Tau Ba 


UX Tau Bb 


0.14 


0.81 


0.57 


<3 


1, 1 


A 




UX Tau Bab 


UX Tau C 


6.45 


0.17 


1.03 


<3 


1, 1 


A 


FX Tau 


FX Tau A 


FX Tau B 


0.89 


0.47 


0.57 


24 ±3 


1, 1 


A 


DKTau 


DK Tau A 


DK Tau B 


2.34 


0.89 


0.64 


57 ±1 


1, 1 


A 


ZZ Tau 


ZZ Tau A 


ZZ Tau B 


0.042 


0.48 


0.40 


<8 


1, 1 


B 


V927 Tau 


V927 Tau A 


V927 Tau B 


0.27 


0.83 


0.40 


< 10 


1, 1 


B 


XZ/HL Tau 


XZ Tau A 


XZ Tau B 


0.29 


0.66 


0.50 


(14 ±3) 


1, 1 


C" 




HL Tau 


XZ Tau AB 


23.31 


0.99 


0.83 


2390 ± 90 


4, 1 


B 


HKTau 


HK Tau A 


HK Tau B 


2.34 


1.00 


0.57 


130 ±3 


1, 1 


A 


V710 Tau 


V710 Tau A 


V710 Tau B 


3.03 


0.70 


0.57 


125 ±3 


4, 1 


A 




V710 Tau AB 


V710 Tau C 


27.97 


0.41 


0.97 


125 ±3 


4, 1 


A 


V827 Tau 


V827 Tau A 


V827 Tau B 


0.093 


0.63 


0.72 


<6 


1, 1 


B 


V826 Tau 


V826 Tau A 


V826 Tau B 


0.014 


1.00 


0.80 


<7 


9== ,9'= 


,B 


V928 Tau 


V928 Tau A 


V928 Tau B 


0.20 


0.95 


0.60 


<8 


1, 1 


B 




V928 Tau AB 


CFHT-Tau-7 


18.25 


0.10 


1.17 


<8 


4, 1 


B 


GG Tau 


GG Tau Aa 


GG Tau Ab 


0.24 


0.83 


0.72 


(557 ±3) 


1, 1 


A 




GG Tau Ba 


GG Tau Bb 


1.48 


0.31 


0.11 


(<5) 


10, lO'' 


,A 




GG Tau Aab 


GG Tau Bab 


10.38 


0.14 


1.32 


(557 ±3) 


4, 1 


A 


FZ/FY Tau 


FZ Tau 


FY Tau 


17.17 


0.89 


0.72 


29 ±7 


4, 1 


B 


UZ Tau 


UZ Tau Ea 


UZ Tau Eb 


0.0009 


0.30 


1.05 


367 ±3 


11" ,11" 


,A 




UZ Tau Wa 


UZ Tau Wb 


0.37 


0.80 


0.50 


65 ±3 


12, 12 


A 




UZ Tau Eab 


UZ Tau Wab 


3.56 


0.66 


1.36 


432 ±4 


4, 1 


A 


JH 112 


JH 112 Aa 


JH 112 Ab 


1.52 


0.02 


0.82 


15 ±4 


1, 1 


A 




JH 112 Aab 


JH 112 B 


6.56 


0.31 


0.83 


15 ±4 


4, 1 


A 


V807/GH Tau 


V807 Tau A 


V807 Tau Bab 


0.28 


0.62 


0.72 


20 ±3 


13, 13 


B 




V807 Tau Ba 


V807 Tau Bb 


0.037 




< 0.45 




13, 


13 




GH Tau A 


GH Tau B 


0.30 


1.00 


0.50 


15 ±3 


1, 1 


B 




V807 Tau AB 


GH Tau AB 


21.77 


0.53 


1.89 


35 ±4 


4, 1 


B 


GI/GK Tau 


GI Tau 


GKTau 


13.14 


1.00 


0.72 


(15 ±1) 


4, 1 


A 


IS Tau 


IS Tau A 


IS Tau B 


0.22 


0.52 


0.64 


30 ±3 


1, 1 


B 


HNTau 


HN Tau A 


HN Tau B 


3.10 


0.27 


0.82 


34 ±3 


4, 1 


A 


IT Tau 


IT Tau A 


IT Tau B 


2.41 


0.29 


0.94 


27 ±4 


1, 1 


A 


J2-2041 


J2-2041 A 


J2-2041 B 


0.42 


0.41 


0.33 


< 19 


1, 1 


B 


HBC 407 


HBC 407 A 


HBC 407 B 


0.14 


0.44 


2.08 


< 14 


1,1 


B 


FF Tau 


FF Tau A 


FF Tau B 


0.037 


0.44 


0.72 


<4 


1, 1 


B 


HBC 412 


HBC 412 A 


HBC 412 B 


0.70 


1.00 


0.50 


<9 


1, 1 


B 


CoKu Tau/3 


CoKu Tau/3 A 


CoKu Tau/3 B 


2.07 


0.17 


0.57 


<8 


1, 1 


B 



33 



Table 6 — Continued 



System 


Pair 


P 


9 


Mp 


-impair 


Refs 




Primary 


Secondary 


["] 




[Mo] 


[mJy] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


HP Tau 


HP Tau/G2 


HP Tau/G3 AB 


10.09 


0.38 


2.49 


«4) 


4, 1, A 




HP Tau/G2 


HP Tau 


21.30 


0.38 


2.49 


(48 ± 2) 






A 




HP Tau 


HP Tau/G3 AB 


16.70 


1.00 


0.95 


(48 ± 2) 






A 




HP Tau/G3 A 


HP Tau/G3 B 


0.030 


0.31 


0.72 


«4) 






A 


Haro 6-28 


Haro 6-28 A 


Haro 6-28 B 


0.65 


0.66 


0.50 


11±3 






B 


HVTau 


HV Tau A 


HV Tau B 


0.036 


0.61 


0.57 


< 7 






A 




HV Tau AB 


HV Tau C 


3.98 


0.78 


0.92 


98 ±3 






A 


VYTau 


VY Tau A 


VY Tau B 


0.66 


0.29 


0.64 


< 10 






B 


GN Tau 


GN Tau A 


GN Tau B 


0.34 


0.87 


0.45 


12 ±3 






B 


JH 223 


JH 223 A 


JH 223 B 


2.07 


0.11 


0.50 


< 7 






B 


IWTau 


IW Tau A 


IW Tau B 


0.29 


0.93 


0.72 


<9 






B 


CoKu Tau/4 


CoKu Tau/4 A 


CoKu Tau/4 B 


0.054 


0.82 


0.54 


9±3 






B 


J04414565+2301580 


J04414565 Aa 


J04414565 Ab 


0.22 


0.10 


0.40 


<24 






A 




J04414565 Aab 


J04414565 B 


12.37 


0.06 


0.44 


<24 






A 


LkHa 332 


LkHa 332/G2 A 


LkHa 332/G2 B 


0.23 


0.64 


0.72 


< 9 




-*- 


B 




LkHa 332/Gl A 


LkHa 332/Gl B 


0.24 


0.55 


0.57 


12 ±3 




-*- 


B 




V955 Tau A 


V955 Tau B 


0.32 


0.63 


0.72 


14 ±2 






B 




LkHa 332/G2 AB 


V955 Tau AB 


10.51 


0.99 


1.18 


14 ±2 






B 




LkHa 332/G2 AB 


LkHa 332/Gl AB 


25.88 


0.75 


1.18 


12 ±3 






B 


DP Tau 


DP Tau A 


DP Tau B 


0.11 


0.74 


0.64 


< 10 






B 


RX J0446.7+2459 


RX J0446 A 


RX J0446 B 


0.051 


0.47 


0.27 


«5) 






D 


DQ Tau 


DQ Tau A 


DQ Tau B 


0.0004 


0.97 


0.65 


208 ±8 


14" , 


14 


B 


Haro 6-37 


Haro 6-37 Aa 


Haro 6-37 Ab 


0.33 


0.11 


0.72 


< 15 






A 




Haro 6-37 Aab 


Haro 6-37 B 


2.62 


0.90 


0.80 


218 ±6 






A 


UY Aur 


UY Aur A 


UY Aur B 


0.88 


0.70 


0.64 


49 ±4 






A 


St 34 


St 34 Aa 


St 34 Ab 


<0.01 


1.00 


0.40 


< 11 


15, 


15 


B 




St 34 Aab 


St 34 B 


1.18 


0.20 


0.80 


< 11 




-*- 


B 


J04554757/J04554801 


J04554757+3028077 


J04554801+3028050 


6.31 


0.70 


0.20 


<23 






A 


HBC 427 


HBC 427 A 


HBC 427 B 


0.032 


0.49 


1.45 


< 14 






B 


CIDA-9 


CIDA-9 A 


CIDA-9 B 


2.34 


0.58 


0.72 


77 ±3 






A 


CIDA-10 


CIDA-10 A 


CIDA-10 B 


0.083 


0.99 


0.27 


< 11 






B 


CIDA-11 


CIDA-11 A 


CIDA-11 B 


0.097 


0.76 


0.27 


< 8 




-*- 


B 


RW Aur 


RW Aur A 


RW Aur B 


1.38 


0.60 


1.20 


71 ±8 






A 



Note. — Col. (1): System name. Col (2): Primary designation of specified pair (note that tliis can consist of more 
than one star; see §2). Col. (3): Secondary designation of specified pair. Col. (4): Angular separation of the pair. 
Col. (5): Mass ratio of the pair [q = Ms/Mp). Col. (6): Estimated stellar mass of the primary component of the pair. 
Col. (7): Composite 880 /im flux density or upper limit for the pair. Values in parenthesis were measured instead at 
1.3mm: the analysis in §5 assumes a conservative scaling by a factor (1.3/0.88)^ ~ 2.2. Col. (8): References for p, q, 
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and Fn; 



Im 



[1] 



Kraus et al. 



2011 



Massarotti et al. (200 5), [10 ] 
[13] = ISchaefer etal] (I2OO6I ) , , 
[B] = [Andrews fc Williamsl (|200 



2] =lBodenet al.l(l2007l ). [3] 



White et al 



14 



(199! 



S 



[11] 



Duchene et al 



Jensen et al 



151 
th| 



(I2OQ: 



2009H), [51 ='Mathieu! (IiQqJ ). [61 =lDudiene et alj (|2006l ) . [71 = ICorrei a et al.' (l2006f )~M 



i), [4] 



(2007), [12] 



Kraus fc Hillenbrand 



Itoh et al.l (1 2005'). [9] = 
Kraus fc Hillenbra nd (20093), 



Mathieu et all (Il997l), [15] =lwhite fc Hillenbrandl (120051). an d [aT^ this article 



[C] = IOsterloh fc Beckwitbl (| 19951 ). [D] = ISchaefer et al.1 (|2009[ ) 



The marginal detections of millimotor emission in the V773 Tau and V410 Tau systems by [Andrews 



Williamg 



(2005J) arc likely un relate d to circumstel lar dust. Both systems are known to emit relatively strong synchrotron radio 
signals iMassi et al.l ( 20061 ) : iBieging et al.l (| 19841 ). To be conservative, we consider these systems to have only upper limits 
on their millimeter luminosities. Synchrotron contami nation from T Tau Sab is also likely to dominate over thermal 
emission from dust, even at millimeter wavelengths (^see iLoinard et al.ll2007l ). The millimeter emission from t he XZ Ta u 



binary may also be contaminat ed, given its relatively bright centimeter- wave signal ICarrasco-Gonzalez et al.l (see 
Moriartv-Schieven et al.l ( 20061 ) suggest a smaller 880 /xm flux density of ~10mJy. 



200' 



i). 



The V773 Tau A-B mass ratio was computed assumi ng t he V773 Tau Aab system mass from lBoden et al.l l|2007f ) 



ng t 
MM 



q. The total mass of L kCa 3 Aab estimated by lKraus et al.l ( 201ll ) was used to help determine the mass r atio of the spec- 
trosco pic binary. The iDuchene et al.l ( 20061 ) results were used in conjunction with the measurements of i White fc Ghed 
(|200ll ) to estimate th e T Tau Sab-N mas s ratio. The DH Tau AB and GG Tau Bab mass ratios implicitly assume the 



Mp values derived bv iKraus et al.l ( 2011). The mass r atio for J 1-4872 Bab was estimated from the spectral type and 



K-band magnitude difference from lCorreia et al.l ( 20061 ). as described bv lKraus et al.l ( 201ll ) 



'^The single-lined spectroscopic binaries LkCa 3 Aab and UZ Tau Eab have a separation weighted by an unknown factor 
of sini and a mass ratio proportional to (si nQ^ (although an independent inclination estimate is available for UZ Tau 
Eab, assuming its gas disk is co-planar; see lSimon et al.l 120001 ). The double-lined spectroscopic binaries V826 Tau AB 



and DQ Tau have only their separations weighted by sini; the mass ratios are determined precisely. IWhite fc Hillenbrand 



((2005) demonstrated that St 34 Aab is a single-lined spectroscopic binary, but did not have sufficient radial velocity 
data to constrain the orbit. Th e V773 Tau B-C proj ected separation was estimated using the weighted mean position 
offsets at the same epochs from lDuchene et al.l (|2003l ). 



Although the current projected separation for the T Tau Sab-N pair is ~0'.'7 (IDuchene et al. 



2OO1 



^, 



astrometric 



monitoring of the pair suggests that the true separation is likely substantially larger. iKohler et al.l ([2008) suggest that 
the pair may be near periastron on an eccentric orbit with a much larger true separation, a ~ 1500 AU, than would be 
suggested by the projected distance on the sky. In the analysis of §5, we adopt their estimate, but note that the major 
features of our statistical analysis do not depend on that decision. 
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Table 7. Single Stars Observed at Millimeter Wavelengths 



Name 


M, 


Fd 


Ref 


Name 


M, 


Fd 


Ref 




[M0I 


[mJy] 






[Mo] 


[mJy] 




(1) 


(2) 


(3) 


(4) 


(1) 


(2) 


(3) 


(4) 


LkCa 1 


0.27 


<8 


1, A 


L1551-55 


0.72 


< 5 


1, A 


IRAS 04108+2910 


0.30 


(<20) 




B 


V830 Tau 


0.72 


<6 




A 


FM Tau 


0.64 


32 ±8 




A 


IRAS 04301+2608 


0.64 


18 + 6 




A 


CWTau 


0.94 


66 + 6 




A 


DL Tau 


0.72 


440 + 40 




A 


FP Tau 


0.27 


«9) 




B 


DM Tau 


0.57 


210 + 1 




E 


CXTau 


0.45 


25 + 6 




A 


CI Tau 


0.72 


324 + 6 




A 


LkCa4 


0.72 


<4 




A 


JH 108 


0.57 


<7 




A 


CYTau 


0.57 


140 + 5 




A 


AA Tau 


0.72 


144 + 5 




A 


V410 Tau X-ray 1 


0.30 


<9 




C 


HO Tau 


0.60 


44 + 6 




A 


HBC 372 


0.94 


<8 




A 


DN Tau 


0.64 


201 + 7 




A 


HBC 376 


0.72 


<6 




A 


HQTau 


0.94 


11 + 3 




A 


BP Tau 


0.72 


130 + 7 




A 


LkCa 14 


0.64 


<9 




A 


V819 Tau 


0.72 


<9 




A 


DO Tau 


0.64 


258 + 42 




A 


DE Tau 


0.50 


90 + 7 




A 


LkCa 15 


1.01 


410 + 1 




E 


RYTau 


1.46 


560 + 30 




A 


ITG 33A 


0.40 


(<5) 




B 


HD 283572 


2.26 


<9 




A 


IRAS 04385+2550 


0.64 


(24 + 1 




B 


FT Tau 


0.60 


121 + 5 




A 


GO Tau 


0.64 


173 + 7 




A 


IRAS 04216+2603 


0.64 


(<20) 




B 


DRTau 


0.72 


533 + 7 




A 


IP Tau 


0.64 


34 + 5 




A 


DS Tau 


0.82 


39 + 4 




A 


DG Tau 


0.94 


1100 + 100 




A 


Haro 6-39 




36 + 6 




, A 


IRAS 04260+2642 


0.94 


(105 + 10) 




D 


CM Aur 


0.94 


640 + 4 




E 


HBC 388 


1.46 


<6 




A 


LkCa 19 


1.69 


< 10 




A 


IQ Tau 


0.60 


178 + 3 




A 


SU Aur 


2.34 


74 + 3 




A 


JH 56 


0.60 


<8 




A 


V836 Tau 


0.72 


74 + 3 




A 


L1551-51 


0.72 


< 13 




A 


CIDA-8 


0.33 


27 + 3 




A 


Haro 6-13 


0.64 


395 + 56 




A 


CIDA-12 


0.27 


<7 




A 



Note. — Col. (1): Star name. Col. (2): Estimated stellar mass. Col. (3): 880 ^m flux density or upper 
limit. Values in parenthesis were measured instead at 1.3 mm: the analysis in §5 ass umes a cons e rvativ e 
scaling by a factor (1.3/0.88)^ Ri 2.2 Col. ( 4): Referen ces for A/, and F^: ill = iKraus et al. (|201ll l, 
and [Al = Andrews fc Williamsl (|2005l ). [B] =ISchaefer et al. (2009.). [C] = lMotte fc Andrei l|200l[| , [D] = 



Mannings fc Emerson 



nl (|l994l ') 



and [E] = [Andrews et all (|201l[ ) 
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Table 8. Millimeter Detection Statistics 



Type 


A''det 


iVtot 


/mm 


singles 


32 


52 


62 ±11% 


binaries 


26-34 


96 


27-35 ± 5% 


triples 


11-15 


39 


28-38 ± 10% 


higher order 


13-18 


44 


30-41 ± 10% 


all multiples 


50-67 


179 


28-37 ± 5% 



Note. — Col. (1): Type of stellar system. 
Col. (2): Number of individual stars with de- 
tected millimeter emission. Col. (3): Total 
number of stars of that category. Col. (4) : Mil- 
limeter detection fraction. 



Table 9. Two-Sample Tests Comparing Millimeter Luminosity Distributions 



Group 1 


Group 2 


P (survival) 


P(K-S) 


single stars 

wide pairs (>300AU) 

medium pairs (30-300 AU) 


wide pairs (>300AU) 
medium pairs (30-300 AU) 
smaU pairs (<30AU) 


0.20-0.40 
<0.00006 
<0.00005 


0.37 
0.017 
0.047 


high mass ratios {q > 0.5) 
wide pairs with q > 0.5 
medium pairs with g > 0.5 


low mass ratios (q < 0.5) 
wide pairs with q < 0.5 
medium pairs with q < 0.5 


<0.04 

<0.00001 

<0.08 


0.48 
0.77 
0.10 


single stars 
primaries in pairs 


primaries in pairs 
secondaries in pairs 


0.41-0.98 
<0.008 





Note. — Cols. (1) & (2): The sub-populations being compared. Col. (3): The 
probability that these two groups are drawn from the same parent distributions. The 
listed values correspond to the range or more conservative limit from the standard 
two-sample tests used in survival analysis (Peto- Prentice, Logrank, and Gehan tests). 
Col. (4): Same as Col. (3), but for the special case that includes only those pairs 
with millimeter-wave detections in each sub-population, derived from the two-sided 
Kolmogorov-Smirnov test. 
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